SPECTRAL DISTRIBUTION STUDIES WITH MODIFIED 
KUO-BROWN INTERACTION 
IN THE UPPER HALF OF THE fp SHELL 

S. Choubey and K.Kar 

Saha Institute of Nuclear Physics 
Block - AF, Sector I, Bidhan Nagar 
Calcutta 700 064, India 

J.M.G. Gomez 

Departamento de Fisica Atomica, Molecular y Nuclear 
Facultad de Ciencias Fisicas, Universidad Complutense de Madrid 

E-28040 Madrid, Spain 

V.R. Manfredi 

Dipartimento di Fisica "G. Galilei", Universita di Padova 
Istituto Nazionale di Fisica Nucleare - Sezione di Padova 
via Marzolo 8, 1-35131, Padova, Italy 

Received: 

ABSTRACT 



The spectral distribution method with modified Kuo-Brown interaction is 
extended to the study of the upper half as well as to odd- A nuclei, of the fp 
shell. The calculations show similar success to that obtained for the lower 
half. 
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Microscopic calculations in the fp-shell involving all four valence orbits 
A/2) /5/2) P3/2 and pi/2 have been actively persued in the last few years. 
Among the residual two body interactions used for this purpose modified 
Kuo-Brown (KB3) has shown remarkable success for the lighter fp-shell nu- 
clei in the full shell model diagonalisation calculations [1,2] as well as in the 
Monte-Carlo shell model studies [3] for both spectra as well as transition 
strengths. Spectral distribution theory originally constructed to reproduce 
the global features of level densities and transition strengths [4] is seen also 
to reproduce binding energies, low lying spectra and transition strengths 
equally well. These studies were performed in detail for the sd-shell [4,5] 
and in some cases for the fp-shell nuclei [6] . Recently such calculations were 
carried out with the KB3 interaction in the lower half of the fp-shell [7] with 
reasonable success. Some of the studies are motivated by applications to 
nuclear astrophysics, such as the problem of calculating the electron cap- 
ture/beta decay rates for supernova and presupernova evolution and r and 
s process nucleosynthesis. Therefore one feels the need to extend the calcu- 
lations to the upper half of the fp-shell as many important nuclei for such 
applications have A>60. We should note here that KB3 was constructed 
to improve the spectroscopy of nuclei at the beginning of the fp-shell and 
in the absence of full diagonalisation the success of KB3 for heavier fp-shell 
nuclei is still uncertain. But KB3 correctly modifies the diagonal matrix ele- 
ments of the 2-body residual interaction and as spectral distribution theory 
relies on the lower order moments of the Hamiltonian in the many body shell 
model spaces, which depend more crucially on the diagonal matrix elements 
than on the nondiagonal ones, we consider it interesting to apply the KB3 
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interaction for the spectral distribution studies in the upper half of the fp- 
shell. In this report we study spectra, occupancies and sum rule strengths 
of transition operators, using the KB3 interaction in the upper half of the 
fp-shell. Here we incorporate the corrections coming from the third and the 
fourth moments of the Hamiltonian in terms of the nonzero skewness (71) 
and excess (72). We shall see that compared to the earlier spectral distribu- 
tion results [6], the agreement of the binding energies with the experimental 
values, particularly for nuclei with A>70, is considerably improved. We also 
apply KB3 interaction to the odd-A nuclei in the fp-shell for the first time 
and see that after taking into account corrections from nonzero (71,72) the 
agreement with experimental data is satisfactory. 

Spectral distribution theory gives smoothed fluctuation free forms for the 
density of states as a function of energy, which in large shell model spaces 
asymptotically go towards gaussians. The formal derivation of the result 
uses the central limit theorem (CLT) for the one body Hamiltonian (H(l)) 
and extends that to the two body Hamiltonian (H(2)), defining an ensemble 
of Hamiltonians and averaging the moments < H p (2) >, p=l,2,...over the 
ensemble [8]. For the 2-body Hamiltonian the ensemble averaged results in 
many particle space follow from the dominance of binary correlations, as 
elegantly demonstrated by Mon and French [8]. Spectral distribution also 
gives polynomial expansions for the expectation values of operators in terms 
of energy where only the first two terms in the expansion contribute in the 
CLT limit [4] . Partitioning the shell model space of m valence particles in N 
single particle states according to configuration and isospin ((m,T) spaces 
where m = mi,m 2 , ...mi are the number of particles in the 1 orbits) and the 
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use of gaussian or gaussian modified by Cornish-Fisher expansion around it 
for the density of states, increases the predictability of the method. To find 
the ground state energy E g one inverts the equation 

£ I"" 9 IrH, T {E)dE = do/2 (1) 

J — oo 

m 

where is the gaussian density of states in the (m, T) spaces normalised 
to d(m,T), the dimensionality of the configuration isospin space, and do 
(=(2J+1)) is the degeneracy of the ground state with spin J. Thus one in- 
tegrates the area below the gaussian configuration densities until the area 
becomes equal to half the ground state degeneracy. This energy value is the 
predicted ground state energy and the method is called the Ratcliff proce- 
dure [9]. As done earlier for the sd-shell [10] and the lower half of the fp-shell 
[7], we improve the predictions by incorporating the corrections from nonzero 
skewness (71) and excess (72) by using the Cornish-Fisher expansion which 
gives 

x = V + j(y 2 - 1) + |(y 3 - 3y) - - 5y) (2) 

where y is the normalised energy (y — (E — e)/cr, e is the centroid and a is 
the width) before the correction and x is the value after it. For the centroids 
and widths we use the values in (m,T) spaces. For 71 and 72, ideally one 
should use the values in (m,T) spaces calculating the third and the fourth 
moments of the (l+2)-body Hamiltonian. But as at present the spectral 
distribution method (SDM) codes can calculate them only in the scalar (m) 
spaces. We use a phenomenological correction term for the excess using 
72 = 72(771) + am + bm 2 where the values of 'a ' and 'b' are obtained through 
a best fit. In the upper half of the fp shell, we find that a parametrised 
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dependence of 72 on isospin makes little improvement in the agreement of 
the calculated binding energies with the experimental ones in contrast to 
the lower half. As the 7 X corrections are very small, we keep the 71 (m) 
unchanged. We calculate the binding energies at fixed configurations before 
and after the corrections and give in Table 1 the difference between the 
predicted and the experimental values given by DIFF= (calculated binding 
energy - experimental binding energy) MeV. The binding energies are taken 
to be positive in agreement with the convention used by experimentalists. We 
also list the corresponding values obtained by Haq and Parikh [6] to show that 
by using the KB3 interaction and by incorporating the (71,72) corrections, 
we obtain considerable improvement over the earlier calculations done by 
them, where they utilise scalar isospin moments with excited state correction 
and use the MHW2 interaction. For the odd-A nuclei the best predictions 
are obtained using a=0.006 and b=-0. 00029 and the RMS deviation of the 
calculated values from the experimental values is 1.22 MeV. For the odd-odd 
nuclei the best predictions are obtained using a=0.005 and b=-0. 00028 and 
the RMS deviation of the calculated values from the experimental values is 
1.34 MeV. For the even-even nuclei the best predictions are obtained using 
a=0.007 and b=-0. 00031 and the RMS deviation of the calculated values 
from the experimental values is 1.77 MeV. One observes that correction of 
72 separately for the odd-A, even-even and odd-odd nuclei brings the values 
closer to experimental ones, as this way one is able to take account of the 
ground state pairing effects to some extent. 

There are indications from experimental pick-up/stripping reaction data 
that the orbit I59/2 starts picking up neutrons when the neutron number 
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of the nuclei goes close to 40 [11]. We have constrained our calculations to 
the four fp-shell orbits and this may be one of the reasons for the somewhat 
larger deviations of the SDM values from the experimental binding energies 
for neutrons almost filling the shell. We note the considerable improvement 
one achieves by indroducing the (71, 72) corrections for these nuclei compared 
to the Haq-Parikh values as seen in Table 1. But we feel future studies for 
A=70 and beyond should include excitations of particles to the 1*79/2 orbit. 

To understand the global properties of the KB3 interaction one should 
study the centroids, widths and its correlation coefficient with other typical 
interactions in the fp-shell evaluated in the (m,T) spaces. For the sake of com- 
parison, we choose the MHW2 interaction [12]. We give the typical example 
of m=24 and m=28 with their two extreme isospin values, but the behavior at 
other particle numbers and isospins is very similar. For m=24, the centroids 
for T=0 and T=8 are -243.89 MeV and -193.21 MeV for KB3 and -229.80 
MeV and -186.49 MeV for MHW2 respectively. For m=28 the centroids for 
T=0 and 6 are -304.12 MeV and -274.55 MeV for KB3 and -285.25 MeV and 
-259.98 MeV for MHW2 respectively.As essentially KB3 differs from MHW2 
in the diagonal two-body matrix elements, one understands the differences 
in the centroids observed. The widths for the MHW2 interaction are seen 
to be different from the KB3 interaction by a factor between 0.75 to 0.90 in 
the upper half. For example for m=24 and T=0 the width for KB3 is 14.89 
MeV, whereas for MHW2 it is 13.09 MeV. As the nondiagonal elements for 
the two interactions are almost identical, the correlation coefficient between 
the two interactions, which has the centroids subtracted out, is always close 
to one throughout the shell. Typically for m=28 it varies between 0.992 to 
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0.989. 

Spectral distribution gives a polynomial expansion in energy for the ex- 
pectation values of operators where the terms beyond the first two in the 
expansion are inhibited by CLT. Explicitly for the operator K in the (m,T) 
space, the CLT form for its expectation value at energy E is, 

K(E; m, T) = (m, T\K\m, T) + <Zk-h(™, T)a K (m, T) E ~ f c ^™\ T > (3) 

where a K (m, T) is the width of the operator K in the space (m,T) and (h-k 
is the correlation coefficient between K and H in the space (m,T). -E c (m,T) 
and a(m, T) are the centroid and width of the Hamiltonian in the same space. 
Taking K — n s , the number operator for the orbit s, one can find the occu- 
pancy of the orbit s. One can also obtain such expansions for the expectation 
values in specific configuration isospin (m,T) spaces and then average over 
all configurations. Table 2 shows such configuration averaged occupancies 
calculated for the 4 orbits at the ground state energies for some typical val- 
ues of valence nucleons and isospins. These results are also available for all 
valence nucleon numbers and isospins. But we mention here that the SDM 
values given are averaged over all J states. Occupancies through shell model 
calculations or from experimental pick- up/stripping sum rules are for spaces 
with fixed J. Ideally one needs to do a J-projection of the SDM value for 
proper comparisons. A detailed comparison of such occupancies by SDM 
and the shell model in the fp-shell, is, in our opinion, of great interest. 

We also use for K the Gamow- Teller (GT) sum rule operator i.e. K = 
(O 11 x O 11 ) 00 with O 11 being the GT transition operator, a vector in both 
J and T. Then for states with isospin zero, we obtain through Eq.(3) the 
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sum rule strength for the Gamow- Teller transition to all final states. In 
Table 3 we give sum rule strengths for self conjugate nuclei (i.e. with N=Z ) 
with valence particles 22, 24, 32 and 34. The table explicitly shows how the 
correlation of K with H changes the sum rule strengths. As the correlation 
coefficient between K and H is very small, the inclusion of the second term 
brings about a decrease of less than 3%. These sum rule estimates of the GT 
strength are useful for the calculations of the electron capture rates on these 
nuclei during the collapse phase of the supernova or the beta decay rates for 
the presupernova evolution [13]. 

One can also use such sum rule estimates for other one body interaction 
operators also. These calculations can be extented to nuclei with nonzero 
ground state isospin. The success of the SDM in reproducing average ener- 
gies as well as transition strengths gives one a method of evaluating many 
structural properties of nuclei with many valence nucleons in active orbits, 
avoids explicit diagonalisation and is also useful for astrophysical applica- 
tions. 

We thank S. Sarkar for his help and for many useful discussions; V. K. 
B. Kota for supplying us some of the SDM programs. This work has been 
partially supported by the Secific Project "Milano 41" and DGES Project 
No. PB96-0604. 
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Table 1 

The difference, DIFF in MeV, between the experimental and the calcu- 
lated binding energies. Column A, B and C give the difference correspond- 
ing to the value calculated by Ratcliff procedure in (m, T) space, by Rat- 
cliff procedure with (7 1 (m),7 2 (m)) correction and by Ratcliff with (71 (m, T), 
72 (m, T)) correction respectively. The last column gives DIFF correspond- 
ing to earlier SDM predictions with excited state correction using the MHW2 
interaction [6]. 



Nucleus 




DIFF 




DIFF 






with KB3 


with 




A 


B 


C 


MHW2 


61 Zn 


18.3 


-0.6 


-0.8 


4.0 


61 Cu 


18.1 


-0.5 


-0.7 


2.0 


61 Ni 


17.6 


0.2 


0.1 


0.9 


Ql Co 


15.1 


1.1 


0.9 


1.7 


61 Fe 


14.4 


2.5 


2.4 


3.1 


63 Ga 


19.2 


1.1 


-0.6 


2.5 


63 Zn 


19.5 


1.5 


-0.1 


3.0 


e3 Cu 


18.0 


1.2 


-0.4 


0.5 




17.7 


2.1 


0.7 


1.8 


63 Co 


16.6 


5.1 


4.0 


1.4 


65 Ge 


19.9 


3.1 


0.1 


0.6 
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Table 1 (contd.) 



Nucleus DIFF DIFF 

with KB3 with 





A 


B 


G 


MHW2 


65 Ga 


19.6 


2.9 


-0.1 


0.4 


e5 Zn 


17.8 


2.6 


-0.1 


-3.9 


65 Cu 


15.9 


2.1 


-0.4 


-5.2 


65 Ni 


12.6 


1.8 


-0.2 


-8.4 


&7 Ga 


17.7 


3.9 


0.0 


-4.6 


67 Zn 


14.6 


3.0 


-0.3 


-10.2 


67 Cu 


11.9 


2.1 


-0.7 


-10.9 


m Ge 


15.9 


4.8 


0.5 


-9.6 


m Ga 


13.7 


3.7 


-0.1 


-13.8 


m Zn 


11.6 


3.2 


-0.1 


-15.0 


m Cu 


4.6 


-0.4 


-2.3 


22.2 


71 Se 


15.7 


6.3 


1.8 


-12.4 


71 As 


13.9 


5.4 


1.3 


-16.6 


G2 Zn 


20.2 


-0.4 


0.1 


2.1 


62 Ni 


18.2 


0.2 


0.7 


0.1 


64 Ge 


19.9 


1.1 


0.0 


2.0 


m Zn 


20.1 


1.3 


0.2 


-0.3 


64 M 


15.9 


1.2 


0.3 


-3.7 
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Table 1 (contd.) 



Nucleus DIFF DIFF 

with KB3 with 





A 


B 


G 


MHW2 


m Ge 


21.2 


3.3 


0.5 


0.7 


m Zn 


18.2 


2.5 


0.1 


-5.2 


m Ni 


12.1 


1.4 


-0.3 


-8.7 


e8 Ge 


19.3 


4.3 


0.4 


-4.3 


68 Zn 


14.1 


2.7 


-0.3 


-12.8 


72 Rr 


8.6 


-0.7 


-5.0 


-20.2 


72 Se 


15.0 


5.7 


1.4 


-14.5 


u Kr 


5.0 


-1.4 


-4.8 


-25.9 


e2 Cu 


19.9 


1.6 


1.1 


4.9 


62 Co 


15.4 


2.7 


0.8 


1.5 


e4 Ga 


23.3 


3.8 


-0.9 


5.6 


e4 Cu 


18.8 


3.0 


-0.8 


1.1 


G4 Co 


14.1 


3.7 


1.2 


2.4 


m Ga 


22.7 


5.5 


-0.3 


2.2 


m Cu 


15.9 


3.7 


-0.4 


-5.0 


e8 Ga 


17.8 


5.5 


0.0 


-6.3 


e8 Cu 


10.8 


3.0 


-0.5 


-14.1 


74 Br 


14.0 


7.9 


3.8 


-18.3 
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Table 2 



Calculated occupancies by SDM for the fp-shell nuclei in the upper half. 



Atomic 


Number of 






Occupancy 




Number 


valence 


Isospin 


A/2 


h/2 


P3/2 


Pl/2 




particles 












64 


24 





14.59 


2.41 


5.20 


1.81 






4 


14.40 


3.24 


4.58 


1.77 


69 


29 


5 


15.38 


4.60 


6.44 


2.59 






11 


15.15 


6.19 


5.38 


2.27 


74 


34 


1 


15.96 


6.62 


7.83 


3.60 






2 


15.00 


3.00 


4.00 


2.00 
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Table 3 



The sum rule strength for Gamow - Teller (GT) transition for T = 
nuclei in fp shell with KB3 interaction by spectral distribution using Eq. (3) 
which includes terms up to CLT(column B). Column A gives values with 
only the first term of Eq.(3). 



no. of GT sum rule strength 



valence 



particles 



A 



B 



22 



14.89 



14.46 



24 



14.44 



14.03 



32 



9.62 



9.44 



34 



7.67 



7.56 
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